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Abstract 

We review the theory for exactly solving quantum Hamiltonian systems through the 
algebraic Bethe ansatz. We also demonstrate how this theory applies to current studies 
in Bose-Einstein condensation and metallic grains which are of nanoscale size. 
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1 Introduction 



The current realisation of nanotechnology as a viable industry is presenting a wealth of chal- 
lenging problems in theoretical physics. Phenomena such as Bose-Einstein condensation, 
entanglement and decoherence in the context of quantum information, superconducting cor- 
relations in metallic nanograins, soft condensed matter, the quantum Hall effect, nano-optics, 
the Kondo effect and Josephson tunneling phenomena are all emerging to paint a vast canvas 
of interwoven physical theories which provide hope and expectation that the emergence of 
new nanotechnologies will be rapid in the short term future. A significant tool in the evo- 
lution of the theoretical aspects of these studies has been the development and application 
of potent mathematical techniques, which are becoming ever increasingly important as our 
understanding of the complexities of these physical systems matures. 

One approach that has recently been raised to prominence in this regard is that of the 
exact solution of a physical model. The necessity of studying the exact solution has been 
demonstrated through the experimental research on aluminium grains with dimensions at the 
nanoscale level. The work of Ralph, Black and Tinkham (RBT) jl| in 1996 detected the 
presence of superconducting pairing correlations in metallic nanograins which manifest as a 
parity effect in the energy spectrum dependent on whether the number of valence electrons on 
each grain is even or odd. A naive approach to theoretically describe these systems is to apply 
the theory of superconductivity due to Bardeen, Cooper and Schrieffer (BCS) |2]. Indeed, the 
BCS model is the appropriate model for these systems, but the associated mean field treatment 
fails. This is because a mean field theory approximates certain operators in the model by an 
average value. At the nanoscale level, the quantum fluctuations are sufficiently large enough 
that this approximation is invalid. In fact, there had been a long harboured notion that 
superconductivity would break down for systems where the mean single particle energy level 
spacing, which is inversely proportional to the volume, is comparable to the superconducting 
gap, as is in the case of metallic nanograins. This was conjectured by Anderson P| in 1959 
on the basis of the BCS theory, but the experiments of RBT show this to not be the case. 
Consequently, an exact solution is highly desired, a view that has been promoted in jlj. 

The study of exact solutions of quantum mechanical models has its origins in the work 
of Bethe in 1931 on the Heisenberg model 0. The field received a tremendous impetus in 
the 1960s with the work of McGuire [B], Yang [7j, Baxter [S] and Lieb and Wu [H], and has 
prospered ever since. The work of RBT cited above has brought the discipline to a new 
audience, when it was realised that the exact solution of the BCS model had been obtained, 
though largely ignored, by Richardson in 1963 JU]- The reason that Richardson's work 
was overlooked for so long is because the theory that had been proposed by BCS was so 
spectacularly successful that there had never been a need to use an alternative approach. 
Once the results of RBT were communicated however, it was clear that a new viewpoint was 
needed. When the condensed matter physics community became aware of Richardson's work, 
his results were promptly adopted and it was shown that the analysis of the exact solution gave 
agreement with the experiments JT]. A concise yet informative account of the developments 
is given in |12j . 

In this review we will recount the quantum inverse scattering method and the associated 
algebraic Bethe ansatz method for the exact solution of integrable quantum Hamiltonians. We 
then show how this procedure can be applied for the analysis of three models which are the 
focus of many current theoretical studies; a model for two Bose-Einstein condensates coupled 
via Josephson tunneling, a model for atomic-molecular Bose-Einstein condensation and the 
BCS model. In each case we undertake an asymptotic analysis of the solution and demonstrate 
how this can be applied to extract the asymptotic behaviour of certain correlation functions 



at zero temperature through use of the Hellmann-Feynman theorem [T!^ . 



2 Quantum inverse scattering method 

First we will review the basic features of the quantum inverse scattering method ^1 [T3] . 
The theory of exactly solvable quantum systems in this setting relies on the existence of a 
solution R{u) G End{V ^V), where V denotes a vector space, which satisfies the Yang-Baxter 
equation acting on the three-fold tensor product space V ® V ® V 

Ri2{u - v)Rn{u)R23{v) = R23{v)Ri3{u)Ru{u - v). (1) 

Here Rjk{u) denotes the matrix in End(y ®V ®V) acting non-trivially on the j-th and fc-th 
spaces and as the identity on the remaining space. The i?-matrix solution may be viewed as 
the structural constants for the Yang-Baxter algebra which is generated by the monodromy 
matrix T{u) whose entries generate the algebra 

i?i2(n - v)T^{u)T2{v) = T2iv)T,{u)Ruiu - v). (2) 

We note that as a result of the Yang-Baxter is necessarily associative. In component form 
we may write 

J^R'kiu - v)mu)T\{v) = Y.Tl{v)n{u)W^{u - v) 
p,<i p,g 

so the Rfj{u) give the structure constants of the algebra. 

Here, we will only concern ourselves with the su{2) invariant i?-matrix which has the form 

R{u) = -^—{u.I(g)I + r]P) 

U + 1] 



/ 1 \ 

b{u) c{u) 

c{u) b{u) 

\ 1 / 



(3) 



with b{u) = u/{u + rj) and c{u) = r]/{u + rj). Above, P is the permutation operator which 
satisfies 

P{x ®y)=y®x ^ X, y E V. 
In this case the Yang-Baxter algebra has four elements which we express as 

Tiu) - ( ^(!> \ . (4) 



C{u) D{u) 

Next suppose that we have a representation, which we denote tt, of the Yang-Baxter 
algebra. For later convenience we set 

L{u)=7,{T{u)) 

which we refer to as an L-operator. Defining the transfer matrix through 

t{u) = tr TT ((T(m))) = TT {A{u) + D{u)) (5) 

it follows from that the transfer matrices commute for different values of the spectral 
parameters; viz. 

[t{u),t{v)]=Q yu,v. (6) 



There are two significant consequences of (jH))- Tlie first is that t{u) may be diagonahsed 
independently of u, that is the eigenvectors of t{u) do not depend on u. Secondly, taking a 
series expansion 

tiu) = CkU^ 
k 

it follows that 

[cfc, Cj] = V/c, j. 

Thus for any Hamiltonian which is expressible as a function of the the operators Ck only, then 
each Ck corresponds to an operator representing a constant of the motion since it will commute 
with the Hamiltonian. When the number of conserved quantities is equal to the number of 
degrees of freedom of the system, the model is said to be integrable. 

An important property of the Yang-Baxter algebra is that it has a co-multiplication struc- 
ture which allows us to build tensor product representations. In particular, given two L- 
operators LF , acting oji V ® U and V ®W respectively, then L = LF is also an 
L-operator as can be see from 

Ri2{u-v)Li{u)L2{v) = Ri2{u-v)L^{u)LYL^{v)LY{v) 

= LUv)L^{u)R,2{u-v)LY{u)LY{v) 

= L^iv)L'^iu)LYiv)LYiu)Ri2iu-v) 

= L]^{v)V^{v)LX{u)LY{u)Ri2{u-v) 

= L2{v)Li{u)Ri2{u- v). 

Furthermore, if L{u) is an L-operator then so is L{u + a) for any a since the i?-matrix depends 
only on the difference of the spectral parameters. 



2.1 Realisations of the Yang-Baxter algebra 

In order to construct a specific model, we must address the question of determining a real- 
isation of the Yang-Baxter algebra. Here we will present several examples which will all be 
utilised later. The first realisation comes from the i?-matrix itself, since it is apparent from (Q) 
that we can make the identification L{u) = R{u) such that a representation of Q is obtained. 
This is the realisation used in the construction of the Heisenberg model [14, 15j. A second 
realisation is given by L{u) = G (c-number realisation), where G is an arbitrary 2x2 matrix 
whose entries do not depend on u. This follows from the fact that [R{u), G ® G] = 0. 

There is a realisation in terms of canonical boson operators b, with the relations [b, 6^] = 
1 which reads JH] 

where = b'^b. There also exists a realisation in terms of the su{2) Lie algebra with generators 

and nacni, 

^ ' u\ -vS- u + r^S' 

with the commutation relations [S'^, S*^] = iS*"*^, [S'"'",S'^] = 25"^. It is worth noting that in 
the case when the su{2) algebra takes the spin 1/2 representation the resulting L-operator 



is equivalent to that given by the i?-matrix. Another is reahsed in terms of the sm(1, 1) 
generators and [HI HHj , 

with the commutation relations [K^,K^] = ifT^, [K^,K~'] = —2K^. 

Below we will use these realisations to construct a variety of exactly solvable models. First 
however, we will introduce the algebraic Bethe ansatz which provides the exact solution. 



3 Algebraic Bethe ansatz method of solution 

For a given realisation of the Yang-Baxter algebra, the solution to the problem of finding the 
eigenvalues of the transfer matrix (jSJ via the algebraic Bethe ansatz is obtained by utilising 
the commutation relations of the Yang-Baxter algebra. We have from the defining relations 
that (among other relations) 

[A{u),Aiv)] = [D{u),D{v)] = 0, 
[B{u),B{v)] = [C{u),C{v)]=0, 

A(u)C(v) = ^^^^^C(v)A(u) ^C(u)A(v), 

U — V u — V 

D(u)C(v) = ^^^^^^C(v)D(u) + -^C(u)D(v). (10) 
u — V u — V 

A key step in successfully applying the algebraic Bethe ansatz approach is finding a suitable 
pseudovacuum state, |0), which has the properties 

A{u) |0) = a{u) |0) , 

fi(n)|0) = 0, 

C(w)|0) ^ 0, 

D{u)\Q) = d{u)\0) 

where a{u) and d{u) are scalar functions. 

Assuming the existence of such a pseudovacuum state, choose the Bethe state 

M 

\v) = \vi,...,VM) = l[C{Vi)\0). (11) 

i=l 

Note that because [C{u), C{v)] = 0, the ordering is not important in (fTTj). The approach of 
the algebraic Bethe ansatz is to use the relations (fTIHl to determine the action of t{u) on \v). 
The result is 

t{u) \v) = A{u, v) \v) 

N / X M 



Eria{vi) YjVi-Vj + rj X 



U - Vi ^f- Vi - Vj 



I ^ rid[vi) T-T fi - - ^ \ I \ i,^. 

+ > 1 1 l^^i, --Vi-i^u, fi+1, vm) (12) 



U - Vi -V^ Vi - Vj 



where 

M M 

A m, v) = a{u) IT + d{u) TT '-. 13 

U — v.: U — Vi 

i=l i=l 

The above shows that 1^) becomes an eigenstate of the transfer matrix with eigenvalue (|13p 
whenever the Bethe ansatz equations 

^=n^^^^^. (14) 

d[Vi) ^f: Vi -Vj + T] 

are satisfied. Note that in the derivation of the Bethe ansatz equations it is required that 
7^ Vj^i, j. This is a result of the Pauli Principle for Bethe ansatz solvable models as 
developed in [12] for the Bose gas. We will not reproduce the proofs for the present cases, as 
they follow essentially the same argument as |19j . 



3.1 Scalar products of states 

One of the important applications of the above discussion is that there exists a formula due 
to Slavnov [T31 1201 HI] for the scalar product of states obtained via the algebraic Bethe ansatz 
for the i?-matrix (jSJ). The formula reads 

S{v:u) = {0\B{vi)...B{vm)C{u,)...C{um)\0) 
detF{u : v) 



where 



detV{u : v) 



d 



dvi ' ' Mj — Vi 



the parameters {vi} satisfy the Bethe ansatz equations (fT^ . and {uj} are arbitrary. The sig- 
nificance of this result is that it opens the possibility to determine form factors and correlation 
functions for any model which can be derived in this manner. Although we will not go into 
any details here, we wish to point out that explicit results for two of the the models which we 
will discuss subsequently can be found in 



4 A model for two coupled Bose-Einstein condensates 

Experimental realisation of Bose-Einstein condensates in dilute atomic alkali gases has stim- 
ulated a diverse range of theoretical and experimental research activity |^ |23 123 12Z1 12H] • A 
particularly exciting possibility is that a pair of Bose-Einstein condensates (such as a Bose- 
Einstein condensate trapped in a double-well potential) may provide a model tunable system 
in which to observe macroscopic quantum tunneling. Below we will show that a model Hamil- 
tonian for a pair of coupled Bose-Einstein condensates admits an exact solution. The model 
is also realisable in Josephson coupled superconducting metallic nanoparticles [221; which has 
applications in the implementation of solid state quantum computers. 

The canonical Hamiltonian which describes tunneling between two Bose-Einstein conden- 
sates takes the form [20] 

H = ^{Nr-N2)'~^iN,-N2)-^ib[h + blh). (15) 
o Z Z 



where b\, denote the single-particle creation operators in the two wells and Ni = b\bi, N2 = 
b\b2 are the corresponding boson number operators. The total boson number A^^i + N2 is 
conserved and set to the fixed value of A^. The physical meaning of the coupling parameters 
for different realisable systems may be found in [211 • useful to divide the parameter space 
into three regimes; viz. Rabi {K/Sj << N~^), Josephson (A^~^ << K/Sj << N) and Fock 
[N << K/Sj). There is a correspondence between (fT5|l and the motion of a pendulum [26] . 
In the Rabi and Josephson regimes this motion is semiclassical, unlike the case of the Fock 
regime. For both the Fock and Josephson regimes the analogy corresponds to a pendulum with 
fixed length, while in the Rabi regime the length varies. An important problem is to study 
the behaviour in the crossover regimes, which is accessible through the exact solution. The 
exact solvability of ([T5|l which we discuss here follows from the fact that it is mathematically 
equivalent to the discrete self-trapping dimer model studied by Enol'skii et al. [30j, who solved 
the model through the algebraic Bethe ansatz. Below we will describe this construction. 

The co-multiplication behind the Yang-Baxter algebra allows us to choose the following 
representation of the monodromy matrix 

L{u) = L\{u + uj)Ll{u-uj) 

{u + uj + riNi){u~uj + riN2) + blbi {u + cu + riNi)b2 + ri-%\ , . 
{u - uj + r]N2)b\ + r]-^ b^ + )' ^ ^ 

Defining the transfer matrix as before through t{u) = ti {L{u)) we have explicitly in the 
present case 

t{u) = + ur]N + rfNiN2 + 7]u{N2 - Ni) + b^bi + b\b2 + r]~'^ - . 



Then 



r]N 



u=0 

and it is easy to verify that the Hamiltonian is related with the transfer matrix t{u) by 
H = -ti (^t{u) - i(t'(0))2 - nt'(O) - r/-2 + 00^ - 

where the following identification has been made for the coupling constants 

K Krf Ajj, Sj 
4 2 ' 2 ' ' 2 

An explicit representation of is obtained from ([T6| with the identification 

A{u) = {u + uj + r]Ni){u-uj + r]N2) + blbi 
B{u) = {u + uj + riNi)b2 + 1]-% 
C{u) = {u - u + r]N2)b\ + r]-^bl 
D{u) = b\b2 + r^. 

Choosing the Fock vacuum as the pseudovacuum, which satisfies B{u)\Q) =0 as required by 
the Bethe ansatz procedure, the eigenvalues a{u) and d{u) of A{u) and D{u) on |0) are 



a{u) 
d{u) 



[u + uj){u — oj), 



The Bethe ansatz equations are then exphcitly 



TV 

with the eigenstates of the form (fTT|) with C{u) given as above. From the Bethe ansatz 
equations, we may derive the useful identity 

m m N 

n'«^'?-^)-n n (^^) 

i=l i=l j=m+l •' 

which will be used later. 

It is clear that the Bethe states are eigenstates of with eigenvalue A^. As A^ is the total 
number of bosons, we expect A^ + 1 solutions of the Bethe ansatz equations. As mentioned 
earlier, we must exclude any solution in which the roots of the Bethe ansatz equations are not 
distinct. For example, the solution 



±^uj'^ - (-1)^77-2, Vj (19) 



of (fT7j) is invalid, except when A^ = 1. (Note the error in [23] )■ For a given valid solution 
of the Bethe ansatz equations, the energy of the Hamiltonian is obtained from the transfer 
matrix eigenvalues (fT^ and reads 

( ^ n rP'N'^ 
E = -K\ri'^T](l + — —)-- ur]N-u'^ 



N 

-r^+^^+(^^-^^)n(i-^) ■ (20) 



Note that this expression is independent of the spectral parameter u which can be chosen 
arbitrarily. The formula simplifies considerably with the choice u = u, hj employing (jlSp . 
which yields a polynomial form. 



^ .2 AT2 



T— I- f] N 

E = — K I r]^^ JJ^ ?72(t'j — + r7)(f j + C(j) rjuN — t]'"^ 

1=1 



However, for the purpose of an asymptotic analysis in the Rabi regime, it is more convenient 
to choose u = 0, while for the Fock regime we use u = rf. 

4.1 Asymptotic analysis of the solution 

Here we will recall the asymptotic analysis of the exact solution that was conducted in (31j . 
We start the analysis with the Rabi regime where TfN << 1. From the Bethe ansatz equations 
it is clear that rfvf — > 1 as — > 0, so that Vi ~ ±77^^. However, when r/ = we know that 
the Hamiltonian is diagonalisable by using the Bogoliubov transformation, from which we 
can deduce that the solution of the Bethe ansatz equations corresponding to the ground state 
must have Vi ~ r]~^. Therefore it is reasonable to consider the asymptotic expansion 

Vi ^ T]-^ + ei + r]6i. (21) 



Excitations correspond to changing the signs of the leading terms in the Bethe ansatz roots. 
To study the asymptotic behaviour for the mth excited state, we set 



Vi ^ -f] ^ + ei + r]Si, i = l,...,m, 

Vi ^ r]-^ + ti + 7]6i, i = m + l,...,N, (22) 



with the convention that the ground state corresponds to m = 0. 

From the leading terms of the Bethe ansatz equations for Vi, i < m we find 



^'-Er^. (23) 

which implies 



m m 



In a similar fashion we have foT m < i < N 



N 



^ 1 



j=7n-\-l 



^ ^ 2 {N ~ m){N - m - I) 



E = 0' E 

i=m+l i=m+l 



which implies 



It is clear from (jzHj) and (|2^ why the Pauli exclusion principle applies in the present case. 
In the asymptotic expansion for fj, is assumed finite. However, if Vi = vj for some i, j, 
then ei = ej and ()23j) and ()24|1 imply that e^, ej are infinite which is a contradiction. Hence Vi 
must be distinct for different i. Note also that for this approximation to be valid we require 
7]^^ >> Cj. However, we see that |ej| is of the order of A^^/^. Thus our approximation will be 
valid for rfN^/'^ << !> which is precisely the criterion for the Rabi region, and consequently 
N cannot be arbitrarily large for fixed r], or vice versa. 

Now we go to the next order. From (fH^ we find 



^ mim — 1) mim — N) muj"^ 



i=l 



^ ^ {N - m){N - m - I) m{m - N) (N - m)uj'^ 
E = i + 2 + 2 

i=m+l 



which using (plj) leads us to the result 



— ^ -N + 2ni- ^ + + —niiN - m). 

K 2 4 2 ^ ' 



The energy level spacings = — E^-i are thus 

f 
2 



^ K(2 + r/V + ^(A^-2m+r 



One may check that Am/N is of the order of A^~^. This indicates that the Rabi regime is 
semiclassical This value for the gap between the ground and first excited state agrees, 
to leading order in ri'^N, with the Gross-Pitaevskii mean-field theory giving a Josephson 
plasma frequency of cuj = 2k(1 + rfN/2Y/'^. 

Now we look at the asymptotic behaviour of the Bethe ansatz equations in the Fock regime 
rf » N . It is necessary to distinguish the following cases: (i) a; = and (ii) u; 7^ 0. 

(i) uj = {). In this case, it is appropriate to consider the permutation operator P which 
interchanges the labels 1 and 2 in (fTH|) . For = 0, P commutes with the Hamiltonian, and 
any eigenvector of the Hamiltonian is also an eigenvector of P with eigenvalue ±1. Therefore 
the Hilbert space splits into the direct sum of two subspaces corresponding to the symmetric 
and antisymmetric wavef unctions. From now on we restrict ourselves to the case when N 
is even, i.e., = 2M, although a similar calculation is also applicable to the case when 
is odd. A careful analysis leads us to conclude that the ground state lies in the symmetric 
subspace. The asymptotic form of the roots of the Bethe ansatz equations for the ground 
state takes the "string"-like structure 

~ _(M - j> ± ^^-^£^r/-(2:'--i) + M(M + l)r^5,^, j = 1, ■ ■ ■ , M 

where is a binomial coefficient. For this asymptotic ansatz to be valid, we require that 
any term in the asymptotic expansion should be much smaller than those preceeding. This 
yields 7f » N which coincides with the defining condition for the Fock region. Throughout, 
the Pauli exclusion principle has been taken into account to exclude any possible spurious 
solutions of the Bethe ansatz equations. 

The above structure clearly indicates that in the ground state the A^ bosons fuse into M 
"bound" states and excitations correspond to a breakdown of these bound states. Specifi- 
cally, the first and second excited states correspond to the breakdown of the bound state at 
— (M — l)rj, with the first excited state in the antisymmetric subspace and the second ex- 
cited state in the symmetric subspace. Explicitly, we can write down the spectral parameter 
configurations for the first two excited states 

vi+ ^ -M?7 + ai+?7"^ vi^ ^ -{M - 1)7] + ai_7]~^ , 



with 



M + 1 _M(M + 1) 

O1+ — , 0-1- — 



a2± 



-(M - If ± (M - 1)V13M2 + 10M+ 1 



12 



(M - 1)(M - 2)v/2M(M + 1) 
" 24 ' 



M - J + 1 

a,± = — , ^ao_i ±, J = 3, ■ ■ ■ , M, 

V(j + l)j(j-l)(j-2) ' ' ' ' 



for the (antisymmetric) first excited state and 



(M + 1)(2M + 1) _ M(M + 1) 
•^1+ ~ 7^ 5 '^1- ~ 



a2± 



-(M - 1)2 ± i{M - l)VllM2 + 14Af-l 
12 ' 



(M - 1)(M - 2)J2M{M + 1) 
asi = ±^ , 

M - J + 1 

a-j± = — a,-i ±, j = 3, ■ ■ ■ , M, 

for the (symmetric) second excited state. The breakdown of the bound state at — (M — 
j = 2, ■ ■ ■ , M results in the higher excited states. 

Substituting these results into (j20|) leads us to the asymptotic ground state energy 



Eo ^ -2/tr7"^M(M + 1), 

while for the first and second excited states we have 

2 _n M2 + M - 2 

El ^ KT] — KT] , 

3 

2 _25M2 + 5M + 2 

E2 ^ KT] + KT] . 

o 

In contrast to the Rabi regime, the Fock regime is not semiclassical, as the ratio of the gap 
A and is of finite order when N is large. 

We can perform a similar analysis for odd A^. In this case, the gap between the ground and 
the first excited states is proportional to ht]"'^ instead of ntf. Furthermore, the ground state 
root structure is different in the odd case since not all the bosons can be bound in pairs. This 
indicates there is a strong parity effect in the Fock regime, in contrast to the Rabi regime. 

(ii) 7^ 0. In this case the root structure is somewhat more complicated than for = 0, 
so we will not present the details. We remark however that our calculations show that up 
to order 77"^ the ground state energy eigenvalue takes the same form as in the case = 0. 
Actually, the leading contribution arising from the u term appears only as w^r/"^. This means 
that the results presented below are applicable for all values of uj (or equivalently A/i). 

Although it is difficult to define rigorously [211 1^^]; the relative phase between Bose- 
Einstein condensates is useful in understanding interference experiments j2Zll2HlEl]- Recall 
that in Josephson's original proposal j35j for Cooper pair tunneling through an insulating 
barrier between macroscopic superconductors, the current is a manifestation of the relative 
phase between the wavefunctions of the superconductors. By definition, the relative phase $ 
is conjugate to the relative number of atoms in the two condensates n = Ni — N2. Using the 
Hellmann-Feynman theorem, we find that 

For the ground state in the limit of strong tunneling (i.e., Rabi regime), < An^ >^ N — 
(A/iiV/^j)^. In the case of weak tunneling (i.e., Fock regime), < An^ 2N{N + 2){Sj/ . 
The degree of coherence between the two Bose-Einstein condensates can be discussed in terms 
of 121 

If t 1 ^^0 

a^ — <a,a2 + a,a, >= 

In the strong coupling limit, a ~ 1 — iV~^(A/i)^/(8£^j)^, indicating very close to full coherence 
in the ground state. In the opposite limit, we have a ~ 2(A^ + 2)£j/K « 1, indicating the 
absence of coherence. The above results give the first order corrections to the results presented 
in |27| for the number fluctuations and the coherence factor at zero temperature. 



5 A model for atomic-molecular Bose-Einstein conden- 



sation 

After the experimental realisation of Bose-Einstein condensation in dilute alkali gases, many 
physicists started to consider the possibility of producing a molecular Bose-Einstein con- 
densate from photoassociation and/or the Feshbach resonance of an atomic Bose-Einstein 
condensate of a weakly interacting dilute alkali gas [^Zl IHH] ■ This novel area has attracted 
considerable attention from both experimental and theoretical physicists, and in particular it 
has recently been reported that a Bose-Einstein condensate of rubidium has been achieved 
comprised of a coherent superposition of atomic and molecular states [2011101 • As stressed in 
IH], even in the ideal two-mode limit, mean field theory fails to provide long-term predictions 
due to strong interparticle entanglement near the dynamically unstable molecular mode. The 
numerical results have shown that the large-amplitude atom-molecular coherent oscillations 
are damped by the rapid growth of fluctuations near the unstable point, which contradicts 
the mean field theory predictions. In order to clarify the controversies raised by these investi- 
gations, one can appeal to the exact solution of the two mode model, the derivation of which 
we will now present. 

The two mode Hamiltonian takes the form 

H = -a^a + —{a''a''b + b^aa), (25) 

where and 6^ denote the creation operators for atomic and molecular modes respectively. 
Note that the total atom number operator N = Na + 2Nf, where Na = a^a, Nh = b^b provides 
a good quantum number since [H, N] = 0. 

In order to derive this Hamiltonian through the quantum inverse scattering method, we 
take the following L-operator 

L{u) = GL\u ~5- 7]-^)L^{u) 

with the matrix G given by 



r]-^ 

This gives us the explicit realisation of the Yang-Baxter algebra 

A{u) = -7]'^ {u + r]K'') {u~6-r]'^ + T]Nb) + bK+ 
B{u) = -K- + Nh) - r]-^b{u - r]K') 
C{u) = 7]-^b\u + 7]K') -r]-^K+ 
D{u) = b'^K' +7]''^{u-r]K') 

and 

t(0) = 6K' + b'^K- + bK+ - r]K'Nb. (26) 

Let |0) denote the Fock vacuum state and let \k) denote a lowest weight state of the 
su(l, 1) algebra with weight k] i.e., \k) = k\k). On the product state |\E') = |0) \k) it is 
clear that B{u)[^) = and 

a{u) = —rj^^{u + rjk){u — 5 — rj^^) 
d{u) = 7]^'^ [u — rfk) . 



We can immediately conclude that the eigenvalues of (|^?)|) are given by 



M M 

A(0) = k{5 + r') n ^ - n ^ (27) 

1=1 1=1 



subject to the Bethe ansatz equations 



{vi + i]k) (1 - r]Vi + r]6) Vi - Vj - rj 



n ^" ^ 

J. J. _ 1, . _|_ n ^ 



{vi - r]k) Vi -Vj + 7] 

Realising the s-u(l, 1) algebra in terms of canonical boson operators through 

2 2 4 

we then find that the Hamiltonian (|^^ is related to (j^UI) through 

E = \imn{t{0) -6 /A) 

with uj = Q6. Note that in this case the possible lowest weight states for the su{l, 1) algebra 
are 

|fc = l/4) = |0), |A; = 3/4) = |0). 

Moreover, we have iV = 2M + 2A; - 1/2. 

It is worth mentioning at this point that another realisation of the sm(1, 1) algebra is given 
in terms of two sets of boson operators by 

= a}c\ K = ac, K = 

2 

with J = Na — Nc a central element commuting with the su{l, 1) algebra in this representation. 
Due to the symmetry ^ we may assume J > 0. For this case we define the Hamiltonian 

H = limfi(t(0) -5/2) +/3J 

= aNa + jN^ + n{a^c^b + b^ac) (29) 

with a = 6Q/2 + (3 and 7 = 6Q/2 — f3. This model has a natural interpretation for atomic- 
molecular Bose-Einstein condensation for two distinct atomic species which can bond to form 
a di- atomic molecule. In this case the possible lowest weight states for the sm(1, 1) algebra 
are 

|A; = (m + l)/2) = (at)"|0) 

and J = 2k — 1. A detailed analysis of this model through the exact solution will be given at 
a later date. 

For the exact solution of the Hamiltonian ()25|) it is necessary to take the quasi- classical 
limit ?7 — i> in the Bethe ansatz equations (j2HI)- The resulting Bethe ansatz equations take 
the form 

S-v, + - = 2y^—. (30) 



Also, in this limit the corresponding energy eigenvalue is 



M 

E = uj{M + k-l/4)-Q'^Vi 

i=l 

AI 

= uj{k-l/A) -2kn^-. (31) 

i=l 

The equivalence of the two energy expressions can be deduced from (j3Up . The eigenstates too 
are obtained by this procedure. Consider the following class of states 

M 

\Vi, ...,Vm) = l[c{v,)\^) (32) 

1=1 

where c{v) = {vb^ - a^ay2), |^) = |0) for A; = 1/4 and |^) = a^O) for k = 3/4. In the 
case when the set of parameters {vi} satisfy the Bethe ansatz equations (jHU)) . then (jH^ are 
precisely the eigenstates of the Hamiltonian. 



5.1 Asymptotic analysis of the solution 

In the limit of large |5| we can perform an asymptotic analysis of the Bethe ansatz equations 
to determine the asymptotic form of the energy spectrum. We choose the following ansatz for 
the Bethe roots 

Vi ~ 5 Vi i < m, 

Vi ^ 5 + ei + 5^^ Hi i > m. 

For i > m we obtain from the zero order terms in the Bethe ansatz equations 

^ 1 



J e,; - e,- 



which implies 



From the terms in 5 we find 



M 



E = 



i=m+l 



M 

/i, = 2(fc + m)+2 ^ 



2 



J — m-f-1 



and thus 

M 

Hi = 2{k + m){M -m). 



i=m+l 

Next we look at the Bethe ansatz equations for i < m. The terms in S give 



2k ^ 1 

l-^i . I . l^j [^i 



which imphes 

m 

^^/ij = —2km — m{m — 1). 



i=l 

This gives the energy levels 



Em ~ ujiM +ik-ll4:))-uiM -m)-VL V V /i^ 

i=m+l 1=1 

= ujim + A; - 1/4) H (3m^ - m + 4fcm - 2fcM - 2mM). 

UJ 

The level spacings are 

A = E — E 1 

^ UJ (M + 2-3m-2A;) 

UJ 

from which we conclude that in this limit the model is semi-classical. 

Let £ denote the ground state energy [£ = Eq for Q6 » 0, £ = Em for Q6 « 0) and 
A the gap to the first excited state. Employing the Hellmann-Ferynman theorem we can 
determine the asymptotic form of the following zero temperature correlations 

^ ^' duj' on 

where 9 = — {a^a)h + 6^aa) is the coherence correlator. For large we introduce the rescaled 
variables 

^• = 7v?7^. ^'-TWr^' W = §^. "-^T^- (33) 

We then have for 5* >> 



while for 5* << 

The above shows that the model has scale invariance in the asymptotic limit. The scaling 
properties actually hold for a wide range of values of the scaled detuning parameter 5*, which 
is established through numerical analysis 



5.2 Computing the energy spectrum 

For this model there is a convenient method to determine the energy spectrum without solving 
the Bethe ansatz equations (cf. JH])- This is achieved by introducing the polynomial function 
whose zeros are the roots of the Bethe ansatz equations; viz. 



M 

1=1 



It can be shown from the Bethe ansatz equations that G satisfies the differential equation 

uG" - {u^ -5u- 2k)G' + {Mu - E/Q + 6{k - 1/4))G = (34) 
subject to the initial conditions 

In order to show this we set 

F{u) = uG" - {u^ -5u- 2k)G'. 

As a result of the Bethe ansatz equations it is deduced that F{vi) = 0. Given that F{u) 
is a polynomial of degree (M + 1), we then conclude that F{u) = {au + [3)G{u) for some 
constants a, /3, which are determined by the asymptotic limits m — > and m — > cx). Eq. 
then follows. 

By setting G{u) = YlnQnU^ the recurrence relation 

_E -u{n + k-l/A) n-M-1 

~ Q{n + l){n + 2k) + + l)(n + 2A;)^""' ^^^^ 

is readily obtained. It is clear from this relation that Qn is a polynomial in E of degree n. We 
also know that G is a polynomial function of degree M and so we must have Qm+i = 0. The 
(M + 1) roots of gM+i are precisely the energy levels Em- Moreover, the eigenstates (jH^ are 
expressible as (up to overall normalisation) 

M 



n=l ^ ^ 



The recurrence relation can be solved as follows (cf. |IHI)- Setting 

n 

9n+i = goYl^jyj, 

j=Q 

with 

_E- w(j + A; - 1/4) 

~ nu + i){j + 2k) ' 

and substituting into the recurrence relation we have 

J-M-1 



' ' (j + l)(j + 2/c)- 
This yields ?/j = 1 + Cj^i/yj^i with 



c 



n\j + l){j + 2k){j-M) 



1 



{E - uj{j + k + 3/4))(E - uj{j + k- 1/4)) - 



which means yj can be expressed as a continued fraction. The requirement that G is a 
polynomial function of order M decrees ?/m = 0, in turn implying 

n^M{M + 2k-l) 



Vm-i 



{E - uj{M + k- 1/4))(E -u{n + k- 5/4)) ' 



which is an algebraic equation that determines the allowed energy levels Em- The above 
procedure can easily be employed to determine the energy spectrum numerically, without 
resorting to solving the Bethe ansatz equations. Explicit results can be found in 



6 The BCS Hamiltonian 



The experimental work of Ralph, Black and Tinkham on the discrete energy spectrum in 
small metallic aluminium grains generated interest in understanding the nature of supercon- 
ducting correlations at the nanoscale level. Their results indicate significant parity effects due 
to the number of electrons in the system. For grains with an odd number of electrons, the 
gap in the energy spectrum reduces with the size of the system, in contrast to the case of a 
grain with an even number of electrons, where a gap larger than the single electron energy 
levels persists. In the latter case the gap can be closed by a strong applied magnetic field. 
The conclusion drawn from these results is that pairing interactions are prominent in these 
nanoscale systems. For a grain with an odd number of electrons there will always be at least 
one unpaired electron, so it is not necessary to break a Cooper pair in order to create an 
excited state. For a grain with an even number of electrons, all excited states have a least one 
broken Cooper pair, resulting in a gap in the spectrum. In the presence of a strongly applied 
magnetic field, it is energetically more favourable for a grain with an even number of electrons 
to have broken pairs, and hence in this case there are excitations which show no gap in the 
spectrum. 

The physical properties of a small metallic grain are described by the reduced BCS Hamil- 
tonian ^T] 

c c 

^ = 5Z ^J'^J ^^Y^ cl^cl_Cj.Cj+. (36) 

Above, j = labels a shell of doubly degenerate single particle energy levels with 

energies ej and rij is the fermion number operator for level j. The operators Cj±, ct^ are the 
annihilation and creation operators for the fermions at level j. The labels it refer to time 
reversed states. 

One of the features of the Hamiltonian (jHUj) is the blocking effect. For any unpaired electron 
at level j the action of the pairing interaction is zero since only paired electrons are scattered. 
This means that the Hilbert space can be decoupled into a product of paired and unpaired 
electron states in which the action of the Hamiltonian on the subspace for the unpaired 
electrons is automatically diagonal in the natural basis. In view of the blocking effect, it is 
convenient to introduce hard-core boson operators bj = Cj_Cj+, fo] = c]_,_c]_ which satisfy the 
relations 

ib]r = 0, [6„ bl] = - 26]6,) [6„ b,] = [6], bl] = 

on the subspace excluding single particle states. In this setting the hard-core boson operators 
realise the su{2) algebra in the pseudo-spin reprepresentation, which will be utilised below. 

The original approach of Bardeen, Cooper and Schrieffer to describe the phenomenon of 
superconductivity was to employ a mean field theory using a variational wavefunction for the 
ground state which has an undetermined number of electrons. The expectation value for the 
number operator is then fixed by means of a chemical potential term fi. One of the predictions 
of the BCS theory is that the number of Cooper pairs in the ground state of the system is 
given by the ratio A/d where A is the BCS "bulk gap" and d is the mean level spacing for the 
single electron eigenstates. For nanoscale systems, this ratio is of the order of unity, in seeming 
contradiction with the experimental results discussed above. The explanation for this is that 
the mean-field approach is inappropriate for nanoscale systems due to large superconducting 
fluctuations. 

As an alternative to the BCS mean field approach, one can appeal to the exact solution of 
the Hamiltonian ()36|) derived by Richardson ^Oj and developed by Richardson and Sherman 
It has also been shown by Cambiaggio, Rivas and Saraceno jl^ that (jHHjl is integrable 



in the sense that there exists a set of mutually commutative operators which commute with 
the Hamiltonian. These features have recently been shown to be a consequence of the fact 
that the model can be derived in the context of the quantum inverse scattering method using 
the L-operator (jHJ with a c-number L-operator |221l^; which we will now explicate. 



6.1 A universally integrable system 

In this case we use a c-number realisation G of the L-operator as well as (jH)) to construct the 
transfer matrix 

t{u) = tro {GoLociu - ec)...Loi{u - ei)) (37) 

which is an element of the £-fold tensor algebra of su{2). Above, trg denotes the trace taken 
over the auxiliary space labelled and G = exp{—aria) with a = diag(l, —1). Defining 

Ti = lim 7r^t(u) 

for j = 1, 2, £, we may write in the quasi-classical limit Tj = tj + o{ri), and it follows from 
the commutivity of the transfer matrices that [tj, r^] = 0, Vj, k. Explicitly, these operators 
read 

r, = 2a5| + ^ (38) 

with 9 = S+^S- + S-^S+ + 2S'^ S'. 
We define a Hamiltonian through 

H - 4 + t ^^±r, -l±C, (39) 

j=i i,fc=i i=i 3=1 

= -ii^^^s^'iiisjst (40) 

i=l j,k=l 

where 

C = S+S- + S-S+ + 2{S'f 

is the Casimir invariant for the su{2) algebra. The Hamiltonian is universally integrable since 
it is clear that [H, tj] = 0, Vj irrespective of the realisations of the su{2) algebra in the 
tensor algebra. 

In order to reproduce the Hamiltonian ()36p we realise the su{2) generators through the 
hard-core boson (spin 1/2) representation; viz 

Sl = b„ ST = b], S]=^-{I-n,). (41) 



In this instance one obtains ()3(j|l (with the constant term — ej) where g = 1/a as shown 
by Zhou et al. [221 and von Delft and Poghossian [13] . 

For each index k in the tensor algebra in which the transfer matrix acts, and accordingly 
in (jlOJ, suppose that we represent the su{2) algebra through the irreducible representation 
with spin s^. Thus {5*^, 5*^, S^} act on a {2sk + l)-dimensional space. In employing the 
method of the algebraic Bethe ansatz discussed earlier we find that 

c 

a{u) = exp[—ari) I I 

u — ek 

k=i 

a[u) = exp[ari) I I 

u — 

k=i 



which gives the eigenvalues of the transfer matrix ()H7j) as 

£ , M 

Am = exp{ar])\\ M 

k=l " j=l J 

C M 

+ exp -ar] M M . 

k=l j=l J 

The corresponding Bethe ansatz equations read 

N rr ^i-^k + ysk vi -vj + y 

exp 2ar/ I I = " I I ■ 

k=i ^k- rjSk vi - Vj -r] 

The eigenvalues of the conserved operators (jHH|) are obtained through the appropriate 
terms in the expansion of the transfer matrix eigenvalues in the parameter rj. This yields the 
following result for the eigenvalues Xj of Tj 

V Mi ^ '=1 ^ V 

such that the parameters Vj now satisfy the Bethe ansatz equations 

2a + V^^ = V^-. (43) 

Through (j42p we can now determine the energy eigenvalues of (|4Up . It is useful to note the 
following identities 

M M C 

2«E^v+2EE^ = A'(^^-i) 

j=i j=i k=i ^ ^ 

M C 

Sk 

ai ■ - 



,M+yy^^ = 

j=l k=l ^ 

M C M C C 

EE;r^-EE^ - ^E- 

j=i fc=i J " j=i k=i ^ " k=i 
Employing the above it is deduced that 



Ea, 


c 

i=i 


2aM 






c 


c 


c c 

j=i k^j 


c 

-2Mj2sk 

k=l 


M 

- 2a ^ Vj + M{M - 1 



which, combined with the eigenvalues 2sj{sj + 1) for the Casimir invariants Cj, yields the 
energy eigenvalues 

M 

E = 2j2^r (44) 



From the above expression we see that the quasi-particle excitation energies are given by twice 
the Bethe ansatz roots {vj} of (jlHI). In order to speciahse this result to the case of the BCS 
Hamiltonian (jHUj) it is a matter of setting = 1/2, V/c. Finally, let us remark that in the 
quasi-classical limit the eigenstates assume the form 

1=1 j = l 

The construction given above can also be applied on a more general level. Taking higher 
spin representations of the su{2) algebra produces models of BCS systems which are coupled 
by Josephson tunneling, as described in ^HlEl- One can also employ higher rank Lie algebras, 
such as so(5) |1H] and sm(4) [49j which produce coupled BCS systems which model pairing 
interactions in nuclear systems. For the general case of an arbitrary Lie algebra we refer to 
|5Uj . Finally, let us mention that if one reproduces the above construction with the su{l, 1) 
L-operator Q in place of the su{2) L-operator (jHJ the pairing model for bosonic systems 
introduced by Dukelsky and Schuck [HI] is obtained. 



6.2 Asymptotic analysis of the solution 

In the limit g ^ we can easily determine the ground state energy of ()36|) : it is given by 
filling the Fermi sea. Below we will assume that the number of fermions is even. Thus for 
small > it is appropriate to consider the asymptotic solution 

Vi ^ ei + g5i + g'^fXi, i = l,...,M. 

Substituting this into and equating the different orders in g yields 

which immediately gives us the asymptotic ground state energy 

2 M C , 



M 2 

2 ^ , - Cfc 



j=l j=l k=M+l 



Next we look at the first excited state. In the g = case this corresponds to breaking 
the Cooper pair at level eM and putting single unpaired electrons in the levels e*/ and eM+i- 
Now these two levels become blocked. To solve the equations (jlHjl for this excited state is the 
same as for the ground state except that there are now (M — 1) Cooper pairs and we have to 
exclude the blocked levels. We can therefore write down the energy 



A/-1 2 ^^-1 

2 ^ , - efc 



Ei^eM + eM+i + 2j2^j-9iM-l) + ^Yl J2 ^ 

j=l j=l k=M+2 

The gap is found to be 



2 

^ V^^Am-e, ^^^^ 



As in previous examples, we can calculate some asymptotic correlation functions for zero 
temperature by using the Hellmann-Feynman theorem. In particular, 



which for i < M gives 



while foT i > M we get 



2 ^ (e, -efc)2 

k=M+l ^ * 



We can also determine the asymptotic form of the Penrose- Onsager- Yang Off-Diagonal Long- 
Range Order parameter [221 ES] to be 



7 t ('1". 



1 dEp 



M C 
j=l k=M+l ^ 



^ ^ ^ , "TT^, 1 ^7 ~ 
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